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Abstract 

We predict a new branch of surface Josephson plasma waves (S JPWs) in layered superconductors 
for frequencies higher than the Josephson plasma frequency. In this frequency range, the permit- 
tivity tensor components along and transverse to the layers have different signs, which is usually 
associated with negative refraction. However, for these frequencies, the bulk Josephson plasma 
waves cannot be matched with the incident and reflected waves in the vacuum, and, instead of 
the negative-refractive properties, abnormal surface modes appear within the frequency band ex- 
pected for bulk modes. We also discuss the excitation of high-frequency SJPWs by means of the 
attenuated-total-reflection method. 

PACS numbers: 74.25.N-, 42.25.Bs 
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High-T c layered cuprate superconductors are important candidates for negative-refractive- 
index (NRI) metamaterials (see, e.g., 

HQ). 

Indeed, being uniaxial strongly anisotropic 
materials, they provide different signs of the permittivity tensor components along, e a b, and 
transverse, e c , to the layers in a wide frequency range (see, e.g., providing a possibility 
of NRI. These metamaterials are attracting considerable attention because they have the po- 
tential for subwavelength resolution and aberration-free imaging. Layered superconductors 
are very promising metamaterials because they are relatively straightforward to fabricate 
(compared to double negative metamaterials) and do not require negative permeability. 

Experiments for the c-axis conductivity in layered superconductors prove the use of a 
model in which the superconducting CuC>2 layers are coupled by the intrinsic Josephson 
effect through the layers (4). Thus, a Josephson plasma with anisotropic current capability 
is produced in layered superconductors. Moreover, the physical mechanisms of the cur- 
rents along and across the layers are fundamentally different. The current along the layers 
is similar to the current in bulk superconductors, while the current across the layers is 
Josephson-type. 

The Josephson current along the c-axis couples with the electromagnetic field inside the 

n 

insulating dielectric layers, forming Josephson plasma waves (JPWs) (see the review |5| and 
references therein). Thus, the propagation of electromagnetic waves through the layers is 
favored by the layered structure. The study of these waves is very important because of 
their Terahertz frequency range, which is still hardly reachable for electronic and optical 
devices. 

Like in common plasma waves, JPWs propagate with frequencies above some threshold 
value (the Josephson plasma frequency uj). However, in the frequency range below uj, the 
presence of the sample boundary can produce surface Josephson plasma waves (SJPWs) jfj], 
which are analog to the surface plasmon polaritons in metals |z, 8]. Such waves can propagate 
along the vacuum-layered superconductor interface and damp away from it. 

At frequencies u higher than the Josephson plasma frequency uj, the normal-to-the-layers 
components of both the group velocity and the Poynting vector of propagating JPWs, have 
signs opposite to the sign of the normal component of the wave- vector k s . This corresponds 
to a NRI. However, the condition u > uj is not sufficient for NRI. This NRI effect can 
be observed at the vacuum-layered superconductor boundary only in a relatively narrow 
frequency range, cuj < cu < w[ ac = uj[e/(e — l)] 1 ^ 2 , where e is the interlayer dielectric 
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constant. A similar limitation exists for any insulator-layered superconductor boundary, if 
the dielectric constant e cxt of the insulator is less than e. The above conditions follow from 
the dispersion relation for the JPWs and the natural limitation for the tangential component 
q (q = ksinO < k, k — ou/c) of the wave- vector for a wave incident at an angle 9 from the 
vacuum onto the layered superconductor. A simple analysis shows that the above inequality 
is compatible with the dispersion relation for JPWs only for frequencies uj < u> < w™ c . In 
other words, any wave with frequency u > u;™ incident from the vacuum cannot propagate 
further in a layered superconductor. JPWs with oj > oj™ c can only match evanescent waves 
in the vacuum. As for frequencies wj < w < u>i ac , the NRI can be observed, but only for 
incident angles 9 higher than some critical value 9 crit - 
For frequencies in the interval 

w[ ac < UJ < LU 2 , 0J 2 = wj7, (1) 

we predict the existence of a new branch of surface waves (here 7 = \ c j X a b 3> 1 is the 
current-anisotropy parameter, A c = c/ujE 1 ' 2 and X a f, are the magnetic-field penetration 
depths along and across the layers, respectively). Despite numerous works on this issue, 
SJPWs with frequencies higher than uj were not discussed before. Here we study them and 
prove that the SJPWs spectrum consists of two branches. The low-frequency branch was 
described in detail in Its spectrum u)(q) follows the "vacuum light line", q = ui/c, and 
deviates from it at frequencies close to uj. The new branch of SJPWs predicted here starts 
at the frequency and follows the vacuum light line for u < w 2 . For frequencies of 

the order of 0J2, the dispersion curve oo(q) strongly deviates from the vacuum light line and 
stops at the frequency u = oj 2 when q w 1/A a &. 

Thus, SJPWs do not exist within the frequency gap ujj < u < w™ c . On the other 
hand, as shown here, this is actually the range where the NRI can be observed (for waves 
incident from the vacuum onto the layered-superconductor boundary). Hence, some kind of 
complementarity between the NRI and surface waves is established in this paper. 

Conditions for the observation of NRI. — Consider a layered structure consisting of su- 
perconducting and dielectric layers with thicknesses s and d, respectively (see Fig. [[]). We 
study the transverse-magnetic JPWs propagating with wave- vector k s = (q, 0, x s ) and hav- 
ing the electric, E s = {E*, 0, E%}, and magnetic, H s = {0,H S ,0}, components proportional 
to exp[z(gx + >c s z — cot)]. The coordinate system is shown in Fig. [TJ 
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FIG. 1: (Color online) Geometry for studying waves in layered superconductors. The interface 
z = divides the layered superconductor from an insulator with dielectric constant e ex t- 

The electromagnetic field inside the layered superconductor is determined by the dis- 
tribution of the gauge invariant phase difference ip(x,z,t) of the order parameter between 
neighboring layers. This phase difference can be described by a set of coupled sine-Gordon 
equations (see review jlj]). In the continuum and linear approximation, tp(x,z,t) can be 
excluded from the set of equations for electromagnetic fields, and the electrodynamics of 
layered superconductors can be described in terms of an anisotropic frequency-dependent 
dielectric permittivity with components e c (Q) and e a b(Q) across and along the layers, re- 
spectively 0. In Ref. 3, the effect of spatial dispersion in e c , related to the capacitive 
interlayer coupling, was taken into account. Here we do not consider this effect because it 
is only important for a very narrow frequency range near ooj . 

In the limit s/d <C 1, the equations for e c (fl) and s a b(fl) can be written as 

1 
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Here we introduce the dimensionless parameters Q = uj/ojj, v ah = Ana a -b/ eujj^ 2 , and v c = 
47ra c /euj. The relaxation frequencies and v c are proportional to the averaged quasi- 
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particle conductivities a a b (along the layers) and a c (across the layers), respectively; coj = 
(STreDjc/hs) 1 / 2 is the Josephson plasma frequency. The latter is determined by the critical 
Josephson current density j c , the interlayer dielectric constant s, and the spatial period of 
the layered structure D = s + d ~ d. 

Analyzing the relations for e c (Jl) and e a b(Q), we conclude that their real parts have 
different signs in a wide frequency range: Uj < uj < U2 = 00 j7 (or 1 < Q < Q 2 = 7 3> 1). In 
this frequency range, the z-components of the group velocity and the Poynting vector of the 
bulk JPWs are directed opposite to the ^-component of the wave- vector k s , and this, at first 
sight, corresponds to having a NRI. However, a more careful analysis shows that this can 
only be observed for a much narrower frequency interval. To verify this, one can consider 
the well-known dispersion relation for the normal component x s of the JPW wave- vector, 

(3) 

that follows directly from Maxwell's equations. Obviously, the JPWs can propagate only 
when Re(x^) > 0, where "Re" stands for the real part. Neglecting dissipation, the permit- 
tivity e a b is negative for the frequency region uj < u < oj 2 considered here. Consequently, 
JPWs can propagate if the factor [k 2 — q 2 /e c (Q)} in Eq. Q is negative. Using Eq. (J2J) for 
e c (Q), one can conclude that this factor is negative only when 1 < Q 2 < 1 + q 2 X 2 . For a 
wave incident, at an angle 6, from the insulator with dielectric constant e cxt onto the layered 

1 /2 

superconductor, we have q = {uej^/c) sin 9. Thus, a NRI can be observed for waves with 
incident angles higher than the critical value 6 C n t defined by the equation, 

sin(0 crit ) = ^e c (tt)/£ cxt . (4) 

It is important to note that, due to the negative sign of the permittivity e a b in Eq. ((3]), the 
incident wave penetrates the superconductor for 9 > 9 CI n, and totally reflects from it for 
6 < #crit; contrary to the standard case of waves incident onto the interface dividing two 
usual right-handed media. 

For e ext < e, Eqs. <^j, (j3J, and the inequality sin(6 l cr it) < 1, provide the conditions, 

( £ V /2 

LOj < UJ < Wi = Uj I . (5) 

\£ — £ cx t / 

Thus, the NRI for a layered superconductor bounded by an insulator with e ext < e can 
only be observed in the frequency range uj < u < u\ = Wjfii. This frequency window 
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can be expanded if one uses an insulator with large enough permittivity e ext . Only for 
insulators with very high permittivity e ext > e, the negative refraction can occur in the whole 
frequency interval uj < 00 < oj 2 - The frequency range for the existence of the bulk JPWs 
in superconductors with capacitive interlayer coupling was derived in [lO]. If the constant 
of this coupling tends to zero, the frequency range obtained in [u| coincides with Eq. §5§. 
Below we consider waves in the insulator-layered superconductor system with e cxt < e for 
frequencies 00 > to\. 

Surface Josephson plasma waves above — When uj < u < u>2 and the factor in Eq. (j3J) 
is positive, the ^-component x s of the wave-vector k s becomes imaginary. This means that 
the wave damps into the layered superconductor. On the other hand, for q > ojy/e ext /c 
the wave damps also into the insulator above the layered superconductor. These are the 
characteristic features of the surface waves discussed in this section. 

Consider an interface (the xy-plane) separating an insulator (z > in Fig. [TJ and a 
layered superconductor (z < 0). We now consider a linear surface transverse-magnetic 
monochromatic wave propagating along the x-axis (i.e., proportional to exp[i(qx — ut)]) and 
decaying into both, the insulator and layered superconductor, away from the interface z = 0. 

Performing the standard procedure for searching surface waves (i.e., solving the Maxwell 
equations for the insulator and layered superconductor with proper boundary conditions 
at the interface between them) we obtain the dispersion relation for the surface Josephson 
plasma waves: 

/ \ N 1/2 

£ext — £<M"J 

'4ct - e c (tt)e ab (tt) ^ 



K (n) = n(e c (n) „ 2 ,^ ) > ( 6 ) 



or, neglecting dissipation, 

1 /2 

Here the dimensionless wave- vector is defined as ft = cq/coj€ cxt . Equation ([7j) describes two 
branches of the dispersion curve for the SJPWs (see Fig. [2]). The first branch exists in the 
low frequency range, < u < loj, and it was studied before in [6[. The second (predicted 

1 /2 

here) branch starts from the light line u = cq/sj^ (or Q = k) at u = u\ (point A in Fig. 
then follows this line, deviates from it at u ~ U2 = 7^j, and stops at the point where 
q = 'jco je 1 / 2 /c, to = CO2 (point B in Fig. EJ. 

Thus, there exists a frequency gap, cjj < cu < cui, in the spectrum of the SJPWs. We 
emphasize that the NRI should only exist within this gap. When the permittivity e cxt of the 

6 
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FIG. 2: (Color online) The dispersion curve (ft = uj/ojj versus n = cq/uje^) for SJPWs at the 
vacuum-layered superconductor interface. The values of the parameters are: 7 = 200, e = 16. 
Inset: zoom-in of the spectrum near the point (k = 1, O = 1). Points A and B correspond to the 
beginning and end of the high-frequency branch. The green dashed line is the vacuum light line 
Q = K. 

insulator increases, the point A in Fig. [2] moves towards point B, and the gap in the SJPW 
spectrum increases. When e cxt = £(1 — 1/7 2 ) ~ e , the points A and B coincide, and the 
high-frequency branch in the SJPWs spectrum disappears. 

Note that the Josephson current is small with respect to the displacement current at 
high frequencies, Q ^> 1. In this case, we can omit 1 in the denominator in Eq. (J7|). 
This corresponds to the dispersion relation for a periodic layered structure without coupling 
between superconducting layers. Specifically, the interlayer Josephson coupling is responsible 
for the appearance of a frequency gap, loj < u < in the spectrum of SJPWs. 

Excitation of the SJPWs above uj\. — It is known that the excitation of surface waves is 
accompanied by the so-called Wood anomalies of the reflectivity and transmissivity coeffi- 
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cients (see, e.g., These resonance anomalies can result in the complete suppression of 
the reflectivity by a proper choice of parameters. Here we consider the excitation of high- 
frequency SJPWs (Q ^> 1) by a wave incident from a dielectric prism with permittivity s p 
onto a layered superconductor separated from the prism by a vacuum gap of thickness 5 (the 
so-called "attenuated-total-reflection" method for producing surface waves, see Fig. [3]). 

The suppression of the specular reflectivity |i?| 2 = \H r /H l \ 2 due to the resonant excitation 
of the surface waves can be observed by changing the incident angle at a given frequency or 
by changing the frequency at a given incident angle, as demonstrated in Fig. [4] (a, b). Here 
H l and H r are the magnetic field amplitudes of the incident and reflected waves, respectively. 
Figure H] (c) shows the sharp decrease of the reflectivity in the (8, Q) plane. 

Conclusions. — Here we predict the existence of a new branch of SJPWs in layered su- 
perconductors for the frequency range higher than the Josephson plasma frequency, which 
is a very unusual phenomenon for plasma-like media. It is important that a NRI can only 
be observed for frequencies within the gap in the spectrum of the SJPWs. Thus, some kind 
of complementarity between NRI and surface waves is established in this paper. We have 
also described the excitation of these SJPWs by means of the attenuated-total-reflection 
method. 

We gratefully acknowledge partial support from the NSA, LPS, ARO, and NSF grant 
No. 0726909. 



[1] M. Ricci et al, Appl. Phys. Lett. 87, 034102 (2005); C. Du et al, Phys. Rev. B 74, 113105 
(2006); B. Wood et al, J. Phys.: Condens. Matter 19, 076208 (2007); F. Magnus et al, Nature 
Materials 7, 295 (2008); E. Narimanov, ibid. 7, 273 (2008); M.B. Romanowsky et al, Phys. 
Rev. A 78, 041110 (2008); K.Y. Bliokh et al, Rev. Mod. Phys. 80, 1201 (2008). 

[2] A.L. Rakhmanov et al. larXiv:0907.3564l (2009). 

[3] V.A. Podolskiy and E. E. Narimanov, Phys. Rev. B 71, 201101 (2005); J.B. Pendry, Science 
306, 1353 (2004); A. Alu et al, J. Opt. Soc. Am. B 23, 571 (2006); O.V. Ivanov et al, 
Crystallogr. Rep. 45, 487 (2000). 

[4] R. Kleiner et al, Phys. Rev. Lett. 68, 2394 (1992); R. Kleiner and P. Muller, Phys. Rev. B 
49, 1327 (1994). 

[5] S. Savel'ev et al. [ arXiv:0903.2969t Rep. Prog. Phys. (to be published). 



8 




FIG. 3: (Color online) A dielectric prism with permittivity e p is separated from a layered supercon- 
ductor by a vacuum gap of thickness 5. An electromagnetic wave with incident angle 8, exceeding 

— 1/2 

the limit angle 9t = arcsin(e p ) for total internal reflection, can excite surface waves that sat- 

1/2 

isfy the following resonance condition: uje p siuO = cq. Here k J and k r are the wave- vectors of 
the incident and reflected waves associated with the magnetic field amplitudes H* and H r . The 
resonance excitation of surface waves by the incident wave produces a strong suppression of the 
reflected wave. 

[6] S. Savel'ev et al, Phys. Rev. Lett. 95, 187002 (2005); V. A. Yampol'skii, et al, Phys. Rev. B 

76, 224504 (2007); V. A. Yampol'skii et al, ibid 79, 214501 (2009). 
[7] P. M. Platzman and P. A. Wolff, Waves and Interactions in Solid State Plasmas (Academic, 

London, 1973). 



9 




incident angle, (deg) 

FIG. 4: (Color online) (a) The reflectivity coefficient |i?| 2 versus the incident angle 9 calculated 
numerically for the parameters v ab = 10~ 3 , 7 = 200, e = 16, e p = 4, and Q = 135. These 
parameters correspond to the solid square on the dispersion curve in Fig. [2j The thickness 5 of 
the vacuum gap (see Fig. [3]) is one wavelength, k5 = 2ir. The vertical dashed line at 9 = 30° 
corresponds to the limit angle of the total internal reflection, (b) The reflectivity coefficient |i?| 2 
versus frequency, f2 = uj/luj, for 9 = 30.6°. (c) Color contour plot of the reflectivity coefficient 
|i?| 2 in the plane (9, fi). The dispersion relation for the waves in the dielectric-vacuum-layered 
superconductor system in Fig. [3] is presented by the solid white curve. 
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